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Abstract 

Using the slender inclusion model developed earlier the 
e 1 as to.s ta t i c interaction problem between a penny-shaped crack 
and elastic fibers in an elastic matrix is formulated. For a 
s ingle set and for multiple sets of fibers oriented perpendic- 
ularly to the plane of the crack and distributed symmetrically 
on concentric circles the problem is reduced to a system of 
singular integral equations. Techniques for the regulariza- 
tion and for the numerical solution of the system are outlined 
'For various fiber geometries numerical examples are given and 
distribution of the stress intensity factor along the crack 
border is obtained. Sample results showing the distribution 
of the fiber stress and a measure of the fiber-matrix inter- 
face shear are also included. 


1 . INTRODUCTION 

The general problem of an elastic matrix containing a 
penny-shaped crack and reinforced by filaments of finite 
length was introduced in [1] and [2]. In this problem the 
external loads were applied to the matrix at "infinity" and 
there were severe stress concentrations around the end points 
of the filaments, meaning that the load transfer from the ma- 
trix to the filaments took place mainly around the ends of the 
filaments. On the other hand, if the medium containing the 
crack is a f i ber- rei nf orced composite with fibers extending 
into the loading grips, outside the stress perturbation zone 

*This work was supported by NS F under the Grant GK-42771 X and 
NASA Langley Research Center under the Grant NGR-39-007-01 1 . 
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of the crack the fiber-matrix interface shear will be zero. 

In this case the region of load transfer between the fiber 
and the matrix will be restricted to the neighborhood of the 
crack (approximately four crack diameters along the fiber), 
and the nature of the governing integral equations will be 
somewhat different than that of the filament problem. 

In this paper the problem of a fiber-reinforced elastic 
matrix containing a penny-shaped crack will be formulated and 
solved under certain simplifying assumptions regarding the 
geometry of the medium. The main assumptions are (a) the 
crack is a plane circular (penny-shaped) internal crack; 

(b) the fibers are circular elastic cylinders which are per- 
fectly bonded to the matrix, are oriented perpendi cul ar to 
the plane of the crack, and are symmetrically distributed on 
circles concentric with the crack; (c) the composite system 
is loaded perpendicular to and away from the crack; (d) the 
fiber diameter is relatively small compared to other lateral 
dimensions so that the assumption of local axisymmetry in 
fiber stresses in considering the fiber-matrix displacement 
compatibility and the use of the filament model developed in 
[1] are justified; and (e) the Poisson's ratios of the matrix 
and the fiber are equal so that the reinforcing fiber may be 
replaced by an auxiliary inclusion having the elastic con- 
stants E^ - E and v and by a layer of body forces [1], where 
E,v are the elastic constants of the matrix and E^ is the 
Young's modulus of the fibers. The technique developed in the 
paper is quite straightforward and may easily be applied to 
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any number of fiber systems each distributed symmetr i cal 1 y 
around the crack. As the results indicate, when the number 
of fibers is increased the stress state around the leading 
edge of the crack approaches an axisymmetric one. Thus re- 
formulating the problem as an axisymmetric one, substantial 
simplification in the solution can be achieved. 

2. FORMULATION OF THE PROBLEM 

First we consider a simple set of fibers distributed 
symmetrically around the crack as shown in Figure 1. As in 
121 the unknown function in the problem is the layer of body 
forces Zfz) imbedded on the fiber-matrix interface which can 
be obtained from the displacement or strain continuity con- 
dition along the interface. If Z{z) is the layer of body 
forces acting on the matrix then the cylindrical auxiliary 
fiber of infinite length is subjected to axial surface trac- 
tions of magnitude -Z(z). Assuming a fixed-grip type of load- 
ing at infinity the equilibrium condition of the auxiliary 
fibergives 

TTr o a o — t — - 2ur /“ Z ( t ) d't - 7T r 2 a f (z) =0 (1) 

t 0 z u 

where r is the fiber radius, cr the stress acting on the ma- 
o o 

f 

trix at z = +», a the stress acting on the auxiliary fiber, 
and E and E the moduli of the fiber and the matrix, respect- 
ively. From (1) the strain of the auxiliary fiber may be 
expressed as 

r Cz(t)dt. < 2 > 
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The strain in the matrix e z may be obtained from the 
superposition of the solutions of four problems shown in 
Figure 2. Thus, the integral equation for Z may be obtained 
by writing 

e z - e z + e z + e z + e z = e z , (-«< z <w). (3) 

Similar to the derivation given in [2], the strains 
3 d 

e z ,..,e z may be expressed as (Figures 1 and 2) 


£ z = ct o /E ' 


(4) 


E z C(z)= HTT £ + k c < z ,t > ] Z C t )dt , ( -~<z«») , 

Zr _o E(k)-1 + 1_ 2r o~ p o 


k c (z ’ t} = p t-z 
0 

2r y( t-z) 


t-z 


[ 2 E ( k ) - K ( k ) 3 


(5) 


n 


1 


+ 7Sr (t-z)Z j 2 3/2 

0 i a 2 [d n . + (t-z)^] J/<: 


(l-2y 


3y ( t- z ) ^ > 

d/ + { t-z ) 2 


Y = 1 / ( 3-4v ) , C 1 = 4w(1-v)/[{l+v)(3-4v)] 

k = 2r /p , p 2 = 4r 2 + (t-z) 2 , 

0 0 0 0 

d i = b [ 2 ( 1 -cos 2lr ^ 1 '~ 1) )] !s , (i = 2....N) 


1+v 


e z ( r ' z ) = - — ( 2 / -rr ) 2 a Q [ ( 1 - 2v ) F q ( r , | z | ) + |z|G(r,|z|)], 


(|z|>0). 


( 6 ) 
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£ z d (r,e,z) 


F 

G 


n 

n 


a 


n 


(r,z) * / 

(r,z) = / 
= 2 -2rN 


= . i±22. { 2 / 7T Z a n b n cos(nNe)* 

*J(1-2v) F n { r , | z | ) + |z|G n (r,|z| )] ,( | z | >0) (7) 

0 V d nN (rs)J nN + 3/ 2 ( s ) e ~ SZds >( z>0 » n = 1 » 2 * • *)( 8 ) 

o c °s 3/2 J nN (rs)J nN+3/2 (s)e 7SZ ds,(z>0 } n = 0 f l,2,..)(9) 

r N ( 2n t +1 ) (-D n ' k do) 

k = ° 2n-2k+l ’ 


b = / a / 2it p ( r , 6 ) rdrde , ■ (11) 

Tia 

b = 2 ( nN+1 ) / a / 2tt p ( r , e ) ( — ) cos (nNe)drde , (12) 

n tt a oo a 

p(r,0) = 47 TTT 7 T /" Z(t)dt !L 4-0 -2v + ^-)dt, (13) 

P m p m 

p 2 m * r 2 - 2br cos(8- ^p)+ b 2 + t 2 (14) 

where a is the radius of the crack which, in this paper, is 
assumed to be the length unit, K ( k ) and E(k) are the complete 
elliptic integrals, N is the number of the equally spaced fi- 
bers around the crack, b is the distance of the fibers from 
the crack center (Figure 1), and (r,6,z) are the cylindrical 
coordinates shown in Figure 1, 0=0 corresponding to fiber 1. 

Substituting now from (2) and (4-7) into (3) with r=b 
and 6*0 we find an integral equation of the following form to 
determine the unknown function Z(z): 

/°° dt + /“ k(z,t)Z(t)dt = f (z) , (-oo< z <co) s (15) 

- oo l - z -0 ° 
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where f ( z ) is a known function and the known kernel k(z,t) 

has an additional logarithmic singularity, after determining 

the body force distribution Z(z), it was shown in [1] that the 

quantities of more physical interest, namely the interface 

shear a and the total fiber stress cr f may be obtained from 
s T 

E 

o s (z) * - E _ E Z ( z ) , (-oo< z <oo), (16) 

a f ( z ) - a f (z) + a z m (z), (-oo< z <cc) (17) 

where a f is the stress in the auxiliary fiber and a m in the matrix. 
The matrix stress is obtained from the superposition shown in 
Figure 3 . Thus , 

a z "(z) = a Q + a z b (b,z) + a z C (b,z) + a z d (b,z), . 

( - co < Z <oo ) ( 1 

where 

a z b (b,z) = - a 0 (2/ir) ! ' 2 [F b (b,|z| ) + | z | G q ( b , | z | ) ] , (19) 


1 -2v 


a z C(b ’ z) = 4'Cl^T Z{t) {t ' z)dt{ [r 2 7 (t . z) S! ] 3/i' 


3 ( t- z ) 


2 . _s2,5/2 ‘ 1 - l r 2 


[r Q c + (t-zTJ 

3 ( t- z ) 2 

[d. 2 + ( t-z ) 2 ] 5/2 


N 

+ Z [ 


1 -2v 


i =2 [d. C + ( t-z ri 


2,3/2 


1 ), 


( 20 ) 


w 

a d (b,z) = - (2/tt)^ £ a n b n [F n (b,|z|) + | z | G n ( b , | z | ) ] . ( 21 ) 

0 

Other quantities of physical interest are the crack open- 
ing displacement u z (r,6, + 0) and the stress intensity factor 
along the leading edge of the crack. Referring to [2] 
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and Figure 2 u z may be expressed as 

u z (r,e,+0} = u 2 b (r, + 0) + u z d (r,9,+0) 



00 

+ Z a b cos ( n N 6 ) H ( r , 0 ) ] , (0<r<l, 
n=0 n n n 

0< 0<2tt ) , 


( 22 ) 


where 

H n (r,z) = /" s'* 5 J nN (rs)J nN+3/2 (s)e' zs ds, (z>0). (23) 

From (see [3]) 


H n (r,+0) = 


we obtain 


n N / , 2.3s 

< r (1-i" ) 

' rz r (3/2) 

0,(1 <r<“) 


(0<r<l ) , 


(24) 


2 CO 

u (r,0,+O) = 4( ,l' v ’ (1-r 2 )^ [c + E a b r nN cos(nN6)] , 

o n = 0 


(0<r<l ,0<6<2ir) . (25) 


Using now the definition of the stress intensity factor as 


k, (6) = lim ( r- 1 ) o m (r,6,0) 
1 r-*-l z 


= - lim — - E ? /2( 1 -r) - u m (r,6,+0), (26) 

r+1 2 ( 1 - v ; 3r z 

we find 


k^e) =(2 /tt)[ o q + Z a n b n cos(nNe)], ( 0 < 9 < 2 rr ) (27) 


Note that because of symmetry the forward and anti-plane shear 
components of the stress intensity factor '(k^ and k^) will 
be zero. 
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3. MULTIPLE SETS OF FIBERS 


The analysis developed for a single set of fibers may 

easily be extended to multiple sets distributed with equal 

angular spacings on concentric rings around the penny-shaped 

crack. Consider the geometry given in Figure 3 and again let 

the fibers be oriented perpendi cul ar to the plane of the crack. 

Let M be the number of concentric rings carrying the fibers, 

b • , ( j = l , . . ,M) the radius of the jth ring, N . , ( j = 1 , . . ,M) the 
J J 

number of the fibers on the jth ring*, and r^ and E j , ( j -1 , . . ,M) 
the radius and the Young's modulus of the fibers in the jth 
set. The unknown functions will now be the layers of body 
forces Z - ( z ) , (j=l,..,M) imbedded in the fiber-matrix inter- 
faces and the strain continuity conditions expressed for a 
representative fiber in each set will give the necessary sys- 
tem of integral equations to determine these function. Thus, 
in this case (dropping the subscript z) the condition (3) is 
modified to read 

£^(z) * bj s 0j ,z) “ e d (z) + ek(bj,z) + £ j C ( bj ,0 j , z ) 

+ effb. ,0. ,z) , -»<z<») (28) 

J J J 


Here it is assumed that the numbers N. are either equal or 
multiples of each other and the fibers are distributed in 
such a way that the body forces (or the fiber stresses) in 
each set is the same. Otherwise, any deviation from symmetry 
will cause the fiber stresses in a given set to be different 
and the total number of unknown functions to be much greater 
than the number of sets M. 
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where 0^ is the angular orientation of the representative 
fiber in the jth set. The strain in the auxiliary fiber, 
e . may be obtained from (2) by replacing e , E^, and Z 

3 z i 

-p 

by e. , E-, and Z,, respectively. Referring to Figures 2 
<3 3 3 

and 3, it is clear;that in the superposition the axisymmet- 

a b 

ric strain components in the matrix e and e (which depend 
only on o Q ) will be the same as that found for the single 
set of fibers in the previous section and are given by (4) 
and ( 6 ) . 


Referring again to Figures 3 and 2c, if we let the 
axial strain in the matrix at a location of a representa- 
tive fiber from the jth set (i.e., at r=b,, 9=0-) due to 

J sJ 

all the body forces Z^(z) imbedded in the fiber-matrix 
interfaces on the ith ring be e?. (b.,0.,z), the third 

J 1 J J 

term in (28) becomes 

M 

e J < b rV z) ■ 6 j i (b j’V z) ’ 


( j =1 , • . »M ; - °° < z < 00 ) . 


(29) 


Similarly, 
d 


£ . 


( b . , 0 . , z ) = E e._. (b . ,0 . ,z) , ( j = 1 , . . , M ; -oo< z <co) ( 30 ) 


3 3 


i =1 


3 1 3 3 


1 J 

where e C .. is obtained from (7-14) by replacing e° z , a , b , 

N, r, 6 in (7-9) by e b i , a jn , b jn , Nj , , 0j ; b o> p, b n> N, 

r Q ,Z in (11-13) by b jo> , b jn , N.. , r. , Zj ; b, 6, N in (14) 

by b., 0-0., N., and cos(nN6) in (12) by cos [ n N . ( 6-0 . ) ] . Note 
3 J 3 J J 

that e c - ■ and are dependent on Z^(z) only. 


To evaluate e • the interaction of two different fibers 
3 

at two arbitrary locations has to be established. We first 
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observe that e^. which appears in {29) is given by (5). For 

J J 

this all one needs to do is to replace s 2 c , b, Z, N, p , and 
c 

r in (5) by e.., b., Z., N., p. and r . , respect i vely . Re- 

v J J J J J J J 

ferring now to Figures 3 and 2c and the basic concentrated 
load solution [1,4], the axial strain at the point r = b . , 

J 

0 = 8., z due to N.. line loads 2irr^Z^(t), acting along r = b.. 
6 = 0 im = 0 i + , z = t, (m = 1,..,N., -oo<t<c») 


may be expressed as 


N i 


r wr i " 1 l-?v 

i • ( b . , e . , z ) = rpi f a Z . { t ) ( t-z ) E { 2 ^ 2-^72 

J1 J J EC l 1 m=1 [d? m + ( t- z) 2 ] 3/2 


3y ( t-z ) 


[d im + (t - z) 1 


2 — 572 ^ dt , (i^j , -«<z<«) (31) 


where 


y = 1 / ( 3- 4v ) , C 1 = 4 tt( 1 - v ) / [ ( 1 + v) (3-4v)], 

d i^ = b j ? + b i 2 - 2 b j b i cos ' 9 i + ' e j>- 


If we now substitute (with the necessary change in notation) 
from(2, 4, 6, 29-31) into (28) we obtain a system of singular 
integral equations of the following form for unknown functions 


» • 


,z 




i 


r~ _sl — 

EC 1 t-z 


Z 4 (t) „ M 

dt + / M E k - • ( z ,t)Z . (t)dt = f i (z), 
i =1 J J 


( j = 1 , . . ,M; -oo< z <co) . 


(32) 


where f. and k . . , (i,j = 1,..,M) are known functions. After 
J J 

solving (32) the interface shear stresses, the fiber stresses, 
the crack surface displacement, and the stress intensity factor 
may be obtained from expressions similar to (16), (17), (with 
18-21), (25) and (27). For example, the stress intensity 
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factor may be expressed as 

M °° 

k, (0) =( 2/tt) {a + Z Z a. b. cos [nN. ( 9-6 .)]},( 0< 0.<2ir ) (33) 

I o j _ i n _Q jn J N J J 

where aj n and b ^ n are obtained from (10-13) by replacing a n , 
b n , N. p, b, 0, Z, r Q by a jn , b jn> H y Pj , b . , 9-8., Z y r y 
respectively. 

4. ON THE SOLUTION OF INTEGRAL EQUATIONS 

The integral equations (15) and (32) represent special 
cases of the following more general system: 

03 / . \ 00 

Ag(x) + dt + / w K( x , t ) g ( t) dt 

/ = f ( x ) , ( - 03 < x < co ) , (34) 

A = (a ij .) # B = ( b i j ) , g ( x ) = (g i (x)) # 

K ( x , t ) = ( k . j ( x , t ) ) , f ( x ) = (f. (x) ) , (i,j = 1 , . . , n ) , 

where the matrices A and B are constant and are such that 

A+B and A-B are nonsingular, the known functions f^ and k^ . 

and the unknown functions g .. may be complex, and the kernels 

k.. satisfy the necessary regularity conditions at infinity. 

* J 

The system of singular integral equations (34) may formally 
be regularized as follows: Define the following holomorphic 

f uncti ons : 


G(z) = 

1 r <» 
2iri -oo 

aill dt. z = x + 

t- Z 

iy , 

(35) 

H( 2) = 

G(Z). 

z cS + 



{ ( A+B ) ' 

1 (A-B)G(z), zeS", 


(36) 

G(z) = 

(G j( z)) 

, H(z) - ( H j ( z ) ) , (j 

= 1 , . . ,n) , 


(A-B)' 1 

= C = 

(c.j). (A+l)' 1 = D = 

(d. .) 
v i J 
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In terms of the sectionally holomorphic functions H.(z) (34) 

J 

may be expressed as 

H + (x) - H"(x) = D[f(x) - /” K(x,t)g(t)dt] . (37) 


The solution of (37) is [5] 

H(z) 3 2 ttT ff ( s ) " K(s ,t)g(t)dt] + P(z)» (38) 

where P(z) = (P.(z)) f (j=l,..»n) is an arbitrary (matrix) poly- 

J 

nomial which must be identically zero because of the condition 
that F(z)+0 for |z|+°°. The solution of (34) may then be ex- 
pressed as 

g ( x ) = G + ( x ) - G"(x) = H + ( x ) - C(A+B)H"(x), (39) 

giving 

g ( x ) + /“ M(x,t)g(t)dt = p(x), (-oo< x <oo) (40) 

where 


M(x.t) = (m jk (x,t)) = ( D + C ) K( x , t ) + C d 5 ’ 


D + C 


D-C 


P ( X ) = ( Pj (x)) = ^ f(x) + ^ f 


Ihl 


°° S -x 


ds 


(j,k = 1 , . . , n ) . (41. a, b) 

It should be emphasized that the objective of the regulari- 
zation process is "smoothing" the kernels. Thus the system of 
singular integral equations (34) is reduced to that of Fredholm 
integral equations given by (40). It is seen that if (34) con- 
sists of its dominant part only, i.e., if K(x,t) = 0, then (40) 
and (41. b) give the solution in closed form. If K is not zero 
and the integral in (40) is one of the standard convolution 



types (see, for example, [6]) the solution of (40) may again 
be obtained in closed form. Otherwise (40) may have to be 
solved numerically. 

From the view point of numerical analysis the procedure 

outlined here which leads to (40) would be very appropriate 

and convenient provided the kernels m . . can be evaluated in 

' 3 

closed form. If these kernels too have to be evaluated numeri- 
cally from the singular integrals given by (41. a), the technique 
could be quite laborious. In this case the following simpler 
and more direct approach may be preferable: Noting that 

iZ 1^7 * 0 * ( -cc<x<co) (42) 

equation (34) may be expressed as 

AgU) + 4- /“ 9( t.x 9(x) dt + C K(x,t)g(t)dt 

= f ( x ) , ( - 00 < x < 00 ) ( 43 ) 

It should be pointed out that because of the nature of the re- 
lated physical problems, unlike the solution of the singular 
integral equations defined on non- i nter secti ng smooth arcs, 
the solution of singular integral equations defined on infinite 
lines and smooth closed contours are usually bounded and con- 
tinuous functions. Hence, for the purpose of numerical analysis, 
in the latter type equations the singularity of the kernel can 
always be removed as in (43) and the resulting system can be 
treated as ordinary Fredholm integral equations. Note that in 
(43) for t=x the integrand in the second term becomes the deri- 
vative of g ( t ) at x which is assumed to be bounded. 
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In solving (43) one may still encounter numerical diffi- 
culties because of the fact that the support of the integral 
equations is infinite. To overcome this usually a simple 
change in variables such as the following would be sufficient: 


t = tan x - tan ^|- s ( -«< ( t ,x ) , 

-1 < ( r , s )<1 ) ( 44 ) 

Also, to improve the effectiveness of the numerical solution, 
it is preferable to use a Gaussi an-type integration formula 
rather than one based on dividing the domain into equal sub- 
intervals. In the type of problems under consideration an 
appropriate integration formula would be the following Gauss- 
Legendre formula: 


f j F ( r ) d r 

where 


n 

E 

j = l 


w j F(r j ) 


(45) 


P n (r j ) = w j “ 


9 i 9 »(j l,..,n). 

{ 1 -r . ) l [P ( r . ) ] 1 
v j 1 n v j 1 J 


After making the transformation (44), (43) becomes 

R 1 9 1 ( r ) ” 9i ( 5 ) j 

A 9l (s) + “ 1 

tan j— - tan cos 2 7p- 

+ f\ K 1 (s , r ) g ( r ) — „ = Ms), (-1<x<1) (46) 

1 cos 2 — 

where g-j(s) = g( tanirs/2) , etc. First it may be noted that the 
regularity condition of the problem requires K(x,t) to decay 
sufficiently fast so that the second integral in (46) is con- 
vergent. Hence, the application of (45) to this integral will 
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be s trai gh tf orward . Unlike the Gauss-Chebyshev integration 

' f 

formulas, in the application of Gaus s-Legendre formula to the 
solution of integral equations both variables s- and r^ are 
taken to be the roots of the samepolynominal, i.e., 

P n (s i ) = 0, P n (rj) = 0, ( i » j = l,..,n). (47) 

this means that a r t i=j the second integral in (46) may need 
some care. Observing that 

g-j (r) - g-j (s) _ 2 dg 1 (r) 

7+ ’n"r i tts . 2 irr " tt dr * 

s+r (tan - tan g— jcos 

this difficulty may easily be circumvented by replacing the 
/ 

integrand in the second term for i=j by 
2 g l^ r i + P ~ ^ 9 1 ^ r i - 1 ^ . 

71 r i+l “ r i-l 

Thus, writing (46) at n locations $ i one obtains a system of 
n (matrix), equations for the unknowns g^s..), (i=l,..,n). 

5. N UMERICAL RESULTS 

( * ) 

Figures 4-13 show some of calculated results. . Figure 
4 shows a typical result giving the body force distribution 
Z(z) for a single f i ber , ( i . e ,N = 1 ) obtained from (15). Note 
that because of symmetry with respect to z = 0 plane Z is an 
odd and the fiber stress a ^ is an even function of z. A 
sample result giving the distribution of the stress intensity 
factor along the crack border for various values of N , the num- 
ber of fibers in a single set is shown in Figure 5. As in the 

““ 

Further results and the details of the analysis may be found 
in [7] . 
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case of filament reinforcements, as N increases the 6-depend- 
ence of the stress intensity factor decreases. Hence, for 
large values of N, for the purpose of analyzing the effect of 
the fiber reinforcement on the stress intensity factor, one 
may treat the problem as being ax 1 symme tr i c . This may be done 
by assuming that the matrix is reinforced by a cylindrical 
membrane having the same modules as the fibers and an equi- 
valent cross-sectional area (i.e., Nut ) . Figure 6 shows the 
result of such an axisymmetric analysis where it is assumed 
that N = 1 6 . It is seen that as the (fiber) distance b increases, 
the effect of the reinforcement rapidly decreases, the stress 
intensity factor ratio k-j / ( 2a o /i"/Tr ) approaching the value for 
the penny-shaped crack without rei nf orcement . For smaller 
values of b as the crack border approaches the fibers the prob- 
lem cannot be treated as axisymmetric regardless of the number 
N. However, the membrane reinforcement has a physical meaning. 
In this case when the crack touches the reinforcement the power 
of stress singularity (see equation (26)) becomes less than -^, 
and therefore the stress intensity factor as defined by (26) 
would be expected to go to zero. This trend is clearly seen 
in Fi gure 6 . 

Figure 7 gives some idea about the reinforcement effect 
of two fibers placed in different locations (b^=1.4a , b^l ,5a , 
6-1=0, 02 = 0. 4tt ) . In this and in the subsequent examples the two 
sets of fibers are assumed to have equal radii ( r -| =r 2 = r 0 ) an< ^ 
sti ff nes s^S( =E 2 = E^) . The effect of the number of fibers in 
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each set is shown in Figure 8 where it is assumed that 
b-|-1.4a, b 2 = 1.5a, N 1 =N 2 = N, r-j=r 2 = 0.2a, E f =130E, and 6^=0, 
6 2 = 'n'/4. It is seen that as N increases the 8-dependence as 
well as the magnitude of the stress intensity factor again 
decreases. Typical results giving the distribution of the 
body forces , Z 2 and fiber stresses cr-^ and for the 
two fibers are shown in Figures 9 and 10. E^=15E used in 

these figures approximately correspond to steel - concrete com- 
bination. Note that at distances beyond a few crack radii 
Z-j and Z 2 die out quite rapidly and and a 2f approach the 

value of a Q E^/E corresponding to the fiber stress in un- 
cracked reinforced matrix. 

In the case of multiple set of fibers too for large 
values of the problem may be treated as axi symme tri c by 
replacing the individual sets of fibers by equivalent elastic 
cylindrical membranes. Figures 11-13 show some of the re- 
sults obtained under this assumption for two sets of fibers. 

In these examples the membranes are assumed to be equivalent 
to 8 fibers in each set. Figure 11 shows the effect of the 
fiber radius r i =r 2 = r o ancl fr ’*9 ure ^ shows the effect of the 
modulus ratio E^/E on the stress intensity factor which is 
seen to decrease with both increasing r Q and increasing E^/E. 
Figure 13 shows the effect of the fiber distance b 2 for a 
fixed b-j on the stress intensity factor k-| . The asymptote 
shown by the dashed line is the stress intensity factor corre- 
sponding to a single set of fibers at b^ . 


It is seen that 



k-j rapidly approaches this asymptotic value as b 2 increases. 
This implies that the sets of fibers which may exist' in the 
composite beyond a certain distance from the crack would have 
only a negligible effect on the stress intensity factor. The 
figure also shows that as b 2 decreases the reduction in the 
stress intensity factor, k-j becomes more significant and again 
k-j would approach zero as the distance from the crack front 
to the reinforcing membrane goes to zero. 

Comparing the results obtained in this paper for the re- 
inforcing fibers with that found in [2] for finite filaments 
it # may be remarked that there is considerable qualitative sim- 
ilarity regarding the distribution of the stress intensity 
factor. However, in the filament case the interface shear (or 
the body force Z) is concentrated around the filament ends 
and generally gives rise to filament stresses which are higher 
than the stresses found in the fibers. As a result the re- 
inforcement effect (measured by the reduction in the stress 
intensity factor) of the filaments is relatively higher than 
that of the comparably distributed fibers. 
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Figure 2. Superposition for the stresses and displacements 
in the matrix. 



Figure 3. Geometry of the multiple set of fibers; M is the 
number of concentric rings containing the fibers; in each 
ring the fibers are equally spaced wi th <*. = 2tt/N . , N * being 

J J J 

the number of fibers in the jth ring. 



Figure 4. Distribution of body force Z(z) for E^=130E, 
v = 0 . 35 , r =0 ,15a, and ( ): b = 1 . 375a, ( ): 





Figure 5. Effect of the number of fibers N on the stress 
intensity factor distribution along the crack border; 
b= 1 . 375a , r Q =0 . 15a, v=0.35, E f =130E. 




Figure 6. Effect of the fiber distance b or the stress in- 
tensity factor for the axisymmetric case; N=16, r o =0.1a, 

E f - 1 30E , v = 0 . 3 5 . 




Figure 7. Angular variation of the stress intensity factor 
for reinforcement by two fibers; b'-j =1 . 4a , b 2 = l . 5a , 6-j = 0 , 

0 2 = O.4tt, E 1 = E 2 = 130E, v = 0 . 3 5 , r -| =r 2 = r 0 5 and { ) : 

r =0 . 1 5a , ( ): r = 0 . 20a (_ ): r=0.30a. 



Figure 8. Angular variation of the stress intensity factor 
for two sets of reinforcing fibers; b-|=1.4a, b 2 = 1.5a, B-|=0, 
e^TT/4, E 1 =E 9 = 130E, v = 0 . 35 , rn=r 9 =0.2a, Nt=N 9 = N=8. 








Figure 9. Distribution of body forces for reinforcement 
by two fibers; b-|=1.5a s b 2 = l . 6a , r^=r 2 = 0.4a, 0 -| = 0 , 0 £ - 0 0 2 tt 

E r E 2 =15E - 



0 2 4 6 8 
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Figure 10. Distribution of fiber stress o f for reinforce- 
ment by two f i bers ; =1 . 5a , b 2 =1.6a, r -| = r 2 = ^*^ a ’ 6^ = 0, 

0 2 = O.2tt, E 1 =E 2 = 15E. 
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Figure 11. Effect of fiber radius r Q =r 1 =r 2 on the stress 
intensity factor for reinforcement by two sets of fibers 
obtained from the axisymmetric solution, N^Ng’8, b ^ = 1 -4a 
b 2 = 1.6a, E 1 =E 2 = 1 30E , v = 0.35. 



Figure 12. Effect of modulus ratio E^/E on the stress in- 
tensity factor for reinforcement by two sets of fibers ob- 
tained from the axi symmetri c solution; E-j=E 2 = E f> v = 0.35, 
N 1 =N 2 = 8 s b 1 -1 . 25a , b 2 = 1.5a, r ] =r 2 = 0.1a. 




Figure 13, Effect of fiber distance b 2 on the stress in- 
tensity factor for reinforcement by two sets of fibers ob 
tained from the axi symme tr i c solution; N.|=N 2 = 8» E-j=E 2 = 130 
o = 0.35, b-, =1 . 25a , r n =r 9 = 0.1a. 



